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PACS 74.50.+r - Tunneling phenomena; point contacts, weak links, Josephson effects 
PACS 74. 25. Op - Mixed states, critical fields, and surface sheaths 
PACS 74 . 72 . Bk - Y-based cuprates 

Abstract. - We report spatially resolved tunneling spectroscopic evidence for field-induced mi- 
croscopic orders in a high-Tc superconductor YBa2Cu307-i. The spectral characteristics inside 
vortices reveal a pseudogap (Vco) larger than the superconducting gap (Asc) as well as a subgap 
(A') smaller than Asc, and the spectral weight shifts steadily from Asc to Vco and A' upon 
increasing magnetic field. Additionally, energy-independent conductance modulations at 3.6 and 
7.1 lattice constants along the Cu-0 bonding directions and at 9.5 lattice constants along the 
nodal directions are manifested in the vortex state. These wave-vectors differ fundamentally from 
the strongly dispersive modes due to Bogoliubov quasiparticle scattering interferences and may 
be associated with field-induced microscopic orders of pair-, charge- and spin-density waves. 



Introduction. — In conventional type-II supercon- 
ductors, superconductivity is suppressed inside periodic 
Abrikosov vortices [1], leading to continuous quasiparti- 
cle bound states and a peak of local density of states 
(LDOS) at zero energy [2-4]. In contrast, the effect 
of magnetic field on high-Tc superconductors is much 
more complicated than that on conventional type-II su- 
perconductors. Macroscopically, high-Tc cuprates are ex- 
treme type-II superconductors with strong thermal, dis- 
order and quantum fluctuations [5-8]. Microscopically, 
neutron scattering experiments on hole-doped cuprate 
Lai.84Sro.i6Cu04 reported an effective radius of vortices 
substantially larger than the superconducting coherence 
length ^sc [9]. Scanning tunneling spectroscopic (STS) 
studies of Bi2Sr2CaCu208+x (Bi-2212) found no zero-bias 
conductance peaks inside vortices [10,11]. Further detailed 
spatially resolved STS studies of Bi-2212 in one magnetic 
field H = 5T revealed a field-induced (4ao x 4ao) conduc- 
tance modulation inside each vortex, where oq — 0.385 
nm is the planar lattice constant of Bi-2212 [12]. The 
latter finding has been attributed to the presence of a co- 
existing competing order (CO) such as pair-density waves 
(PDW) [13,14], pinned spin-density waves (SDW) [15,16], 
or charge-density waves (CDW) [17-21] upon suppres- 



sion of SC inside the vortices. However, there have not 
been high-resolution STS studies of the field-induced col- 
lective modes on other high-Tc superconductors than Bi- 
2212, neither has there been detailed investigation of H- 
dependent evolution of the collective modes, although vor- 
tex images had been reported in Y-123 [22,23] and effects 
of quasiparticle scattering by vortices have also been in- 
vestigated in Ca2-xNaxCa02Cl2 [24]. 

In this letter we report spatially resolved vortex-state 
STS studies of Y-123 as functions of magnetic field, which 
reveal two field-induced energy scales and three density- 
wave modes with energy-independent wave- vectors, in ad- 
dition to the effects associated with quasiparticle scatter- 
ing interferences [25,26]. These field- induced microscopic 
orders differ fundamentally from the predictions due to 
simple Bogoliubov quasiparticle scattering [25,26] and are 
suggestive of significant interplay between SC and com- 
peting orders upon increasing magnetic field. 

Experimental. — The synthesis and characterization 
of the Y-123 untwinned single crystal with Tc = 93 K 
studied here have been described elsewhere [27]. We use 
chemical etching to prepare Y-123 surface for STS exper- 
iments because Y-123 is more difficult to handle than Bi- 
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2212 [23,28], and bromine chemical etching techniques can 
reproducibly remove non-stoichiometric surface layers to 
reveal quality surfaces, as manifested by x-ray photoemis- 
sion spectroscopy [28] . 

The spatially resolved tunnelling conductance {dl/dV) 
versus energy (lu = eV) spectra for the quasiparticle 
LDOS maps were obtained with our homemade cryogenic 
scanning tunnelling microscope (STM). Our STM has a 
base temperature of 6 K, variable temperature range up 
to room temperature, magnetic field range up to 7 Tesla, 
and ultra-high vacuum capability down to a base pressure 
< 10~^ Torr at 6 K. For each constant temperature (T) 
and magnetic field (H), the experiments were conducted 
by tunneling currents along the crystalline c-axis under 
a range of bias voltages at a given location. The typical 
junction resistance was ~ 1 Gfi. Current (/) vs. voltage 
{V) measurements were repeated pixel-by-pixel over an 
extended area of the sample. To remove slight variations 
in the tunnel junction resistance from pixel to pixel, the 
differential conductance at each pixel is normalized to its 
high-energy background [29] . 

Results and Analysis. — In Fig. 1(a) we illustrate 
the normalized zero-field c-axis tunneling conductance 
spectra taken at T = 6, 77 and 102 K. For T = 6 K, 
the spectrum exhibits clear coherence peaks at energies 
±Asc ~ ± 20 meV and shoulder-like satellite features at 
zbAoff ~ ± 38 meV. At T = 77 K < T^, only one set 
of rounded features remains. Eventually for T = 102 K 
> Tc the peaks vanish within experimental resolution. The 
long-range homogeneity of the zero-field tunneling spec- 
tra is exemplified by the (95 x 95) nm^ spatial map of 
the superconducting (SC) gap Asc in Fig. 1(b) and by 
the corresponding energy histogram (Asc = 20.0 ±1.0 
meV) in Fig. 1(c), which differs from the strong spatial 
variations in the quasiparticle spectra of Bi-2212 [25]. In 
contrast, the satellite features at ±Aoff exhibit stronger 
spatial inhomogeneity, as manifested by the (95 x 95) nm^ 
spatial map of Aeft in Fig. 1(e) and by the corresponding 
histogram in Fig. 1(f), showing Agg = 37.8 ±1.2 meV. 

A possible interpretation for the zero-field spectra in 
Fig. 1(a) may be a scenario of coexisting CO and SC in 
the ground state of the cuprates [19-21,26,30]. Follow- 
ing the analysis briefly outlined in the Discussion section 
and detailed elsewhere [19,20], we can account for the 
zero- field spectra in Fig. 1(a) by incorporating realistic 
bandstructures for Y-123 and assuming either CDW or 
disorder-pinned SDW as the CO with a wave-vector Qco 
parallel to the Cu-0 bonding directions. We find that 
both CDW and disorder-pinned SDW yield equally good 
fitting except [QcdwI = 2|QsdwI [16]- At T < our 
theoretical fitting (black solid curve) can account for the 
sharp superconducting coherence peaks at lu — ±Asc a-nd 
the shoulder-like satellite features at ui = ±Aeff, where 
the effective gap A^s is related to Age and the CO en- 
ergy Vco via the relation A^^^ = [A^^ + V^^]. [19,20] 
Thus, we obtain Asc = 20 meV, Vco = 32 meV and 
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Fig. 1: (color online) Implication of CO from zero- and finite- 
field STS in Y-123: (a) Normalized zero-field tunneling spectra 
taken at T = 6, 77 and 102 K. The solid lines represent fittings 
to the T = 6, 77 and 102 K spectra by assuming coexisting SC 
and CO, following Refs. [19,20]. For more details of the data 
normalization and theoretical fitting procedures, see Discussion 
and Refs. [19,20]. (b) The Asc map over a (95 x 95) nm^ area 
at T = 6 K and H = 0. (c) Histogram of Asc over the same 
area as in (b), showing Asc = (20.0 ± 1.0) meV. (d) Spatially 
averaged intra- and inter-vortex spectra for T = 6 K and H = 
2 T. (e) The Aefi map over a (95 X 95) nm^ area at T = 6 K 
and H = 0. (f) Histogram of Acff in the same area as in (e), 
showing Acff = 37.8 ± 1.2 meV. 

IQcol — (0.25 ± 0.03)7r. For elevated temperatures be- 
low Tc, such as r = 77 K, our fitting (red solid curve in 
Fig. 1(a)) also agrees with experimental data with rounded 
features at ±Aeff(T) [19,20]. 

An alternative way to verify the feasibility of the CO 
scenario is to introduce vortices because the suppression 
of SC inside vortices may unravel the spectroscopic char- 
acteristics of the remaining CO. As exemplified in Fig. Id 
for a set of intra- and inter-vortex spectra taken at iJ = 2 
T, the quasiparticle spectra near the center of each vortex 
exhibit pseudogap (PG)-like features, in contrast to theo- 
retical predictions for a sharp zero-energy peak around the 
center of the vortex core had superconductivity been the 
sole order in the ground state [2,3,31]. Interestingly, the 
PC energy inside vortices is comparable to the CO energy 
Vco ~ 32 meV derived from our zero- field fitting as well 
as the spin gap energy obtained from neutron scattering 
studies of optimally doped Y-123 [32]. 

To investigate how quasiparticle spectra evolve with 
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Fig. 2: (color online) Vortex-state conductance maps at T = 
6 K in Y-123: (a) Conductance power ratio vq map over a 
(75 X 38) nm^ area for H — 2 T, showing disordered vortices 
with an average as = (33.2 ± 9.0) nm. (b) The re map over a 
(75x40) nm^ area for H = 4.5 T, showing as = (23.5±8.0) nm. 
(c) Conductance spectra along the white line in (a), showing 
SC peaks at uu — ±Asc outside vortices and PG features at 
LJ = ±Vco inside vortices, (d) Conductance spectra along 
the dashed line indicated in (b). Spatially averaged intra- and 
inter-vortex spectra for (e) H ^ 2 T, (f) H ^ 4.5 T and (g) H 
= 6 T. 



field, we performed spatially resolved spectroscopic stud- 
ies at H = 2, 4.5, 5 and 6 T and for T = 6 K. In Figs. 2(a) 
and 2(b) we show exemplified spatial maps of the conduc- 
tance power ratio ra for H = 2 and 4.5 T, respectively. 
Here ra at every pixel is defined as the ratio of the con- 
ductance power [dl /dVY' at ut — Asc relative to that at 
oj = 0. We find that the presence of vortices is associ- 
ated with the local minimum of re because of enhanced 
zero-energy quasiparticle density of states inside the vor- 
tex core. Moreover, the total flux is conserved within the 
area studied despite the appearance of disordered vortices. 
That is, the total number of vortices multiplied by the flux 
quantum is equal to the magnetic induction multiplied by 
the area, within experimental errors. Thus, we obtain av- 
eraged vortex lattice constants = 33.2 nm and 23.5 nm 
for H — 2 T and 4.5 T, respectively, comparable to the 
theoretical values of qb = 35.0 nm and 23.3 nm. On the 
other hand, the mean "vortex halo" radius ^haio is much 
longer than the SC coherence length ^sCi and the average 
Chaio decreases with field. We find ^haio — (7.7 ± 0.3) nm 
ioT H ^2 T, (6.4 ±0.6) nm for H = 4.5 T, and (5.0 ±0.7) 
nm for iJ = 6 T. 

The spatial evolution of the vortex-state spectra may 



be better manifested by following a line through multiple 
vortices in the vortex maps of Figs. 2(a) and 2(b). As 
shown in Figs. 2(c) and 2(d) for H = 2 and 4.5 T, respec- 
tively, the vortex-state spectral characteristics differ from 
those of conventional type-II superconductors [4], show- 
ing modulating gap-like features everywhere without any 
zero-energy peaks. In Figs. 2(e) - 2(g) we compare repre- 
sentative spectra taken inside and outside of vortices for 
H = 2, 4.5 and 6 T. In a constant field, the inter-vortex 
spectrum reveals a sharper set of peaks a.t lu — ±AsCj 
whereas the intra-vortex spectrum exhibits PG features 
at w = ±Vco and Vco > Asc- Additional subgap fea- 
tures at LU — ±A' — ±(7 ^ 10) meV are found inside 
vortices, which become more pronounced with increasing 
H. The physical origin of A' is still unknown. 

Next, we perform Fourier transformation (FT) of the 
vortex-state spectra and compare the results with the 
FT zero-field spectra taken in the same area. We de- 
fine the FT-LDOS taken under a magnetic field H and 
at a constant energy u by the quantity F i\<i, uj , H) . In 
Figs. 3(a) - 3(c), we illustrate field-induced FT-LDOS 
|F(k,w,iJ)| = \F{\l,uj,H) - i^(k,u;,0)| at = 5 T and 
integrated over different ranges of energies, where 



F(k,c^,i/) 
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dV 
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Here denotes the coordinate of the i-th pixel, and 
the sum is taken over all pixels of each two-dimensional 
map. Systematic analysis of the energy dependence of 
the FT-LDOS [33] reveals two types of diffraction spots. 
One type of spots are strongly dispersive and may be at- 
tributed to elastic quasiparticle scattering interferences 
as seen in the zero-field FT-LDOS [25,26]. The other 
type of spots are nearly energy-independent, as mani- 
fested in Figs. 3(a) and 3(b) for FT-LDOS at two dif- 
ferent energies u —10 meV and —20 meV, respec- 
tively, and in Fig. 3(c) for normalized FT-LDOS inte- 
grated from LU = —1 meV and —30 meV. In addition to 
the reciprocal lattice vectors, we find two sets of nearly 
energy-independent wave-vectors along (tt, 0)/(0, tt) and 
one set along (tt, tt): Qpdw = [±(0.56 ± 0.06)7r/ai, 0] and 
[0,±(0.56± 0.06)7r/a2], Qcdw = [±(0.28 ± 0.02)7r/ai, 0] 
and [0,±(0.28 ± 0.02)7r/a2], and Qg^w = [±(0.15 ± 
0.01)7r/ai,±(0.15 ± 0.01)7r/a2]. Here ai = 0.383 nm 
and 02 = 0.388 nm. For clarity, we illustrate |-F'(k, w)! 
for different energies along k || (7r,0) in the upper of 
Fig. 3(d) and along {tt,tt) in the lower panel, where the 
w-independent peaks correspond to [QxdwI (^ = P; C, 
S) and the reciprocal lattice vector (27r/ai_2)- These 
wave- vectors are nearly energy independent, as shown in 
Fig. 3(e). For comparison, a dispersive wavevector due to 
quasiparticle scattering interferences [25,26] along (tt, tt), 
which is denoted as 97 in Ref. [25], is also shown in the 
lower panel of Fig. 3(e). We note that the dispersion rela- 
tion for the mode along the nodal direction is in good 
agreement with both the experimental results found in Bi- 
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Fig. 3: (color online) FT studies of the vortex-state conductance maps in the two-dimensional reciprocal space and at _ff = 5T: 
(a) FT-LDOS \F{\^,uJ)\ for oj = -10 meV. (b) FT-LDOS \F{^^,uJ)\ for uj = -20 meV. (c) Normalized FT-LDOS obtained by 
integrating [^(k, t<j)| from —1 meV to —30 meV. Comparing (a) - (c), we find three sets of energy- independent spots in addition 
to the reciprocal lattice constants and the (tt, tt) resonance: Qpdw ^'^'^ Qcdw along the (tt, 0) / (0, tt) directions and Qsdw along 
(tTjTt), which are circled for clarity, (d) |-F(k,aj)| for different energies are plotted against k || (7r,0) and (vr,7r) in the upper and 
lower panels, respectively, showing peaks at energy-independent Qpdwi Qcdw and the reciprocal lattice constants at (27r/ai) 
along (7r,0), and at Qsdw along {n,-n). Additionally, dispersive wavevectors due to quasiparticle scattering interferences are 
found, as exemplified in the lower panel, (e) Momentum (|q|) vs. energy (uj) for [Qpowli IQcdwI and IQsdwI- O^i® dispersive 
wavevector along [11,11), denoted as qj [25], is also shown in the lower panel for comparison, (f) The symmetric and anti- 
symmetric components of Re[_F(k, ci;)] for k = Qpdwi Qcdw and Qsdw are shown as functions of to in the upper and lower 
panels, respectively. 
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2212 [25] and the theoretical predictions for quasiparticle 
scattering interferences [25,26]. 

Interestingly, we find that Qp^w and Qcdw correspond 
to charge modulations at wavelengths of (3.6±0.4)ai^2 and 
(7.1 ± 0.6)ai.2 along the Cu-0 bonding directions. The 
former is comparable to the checkerboard modulations re- 
ported in the vortex state of Bi-2212 [12], whereas the lat- 
ter is consistent with the Fermi surface-nested CDW wave- 
vector derived from our zero-field analysis. On the other 
hand, the nodal wave-vectors Qsdw niay be associated 
with SDW because of the field-induced unequal popula- 
tions of spin- up and spin-down quasiparticles. In this con- 
text, we note that intense spots associated with the (tt, tt) 
spin resonance [32] are also manifested for ui < A^c, as ex- 
emplified in Fig. 3(a) for cj = —9 meV. These field- induced 
collective modes along the nodal direction are suggestive 
of important interplay between SC and spin excitations in 
the cuprates. 

To better understand the nature of these field-induced 
wave- vectors, we consider the symmetry of the complex 
quantity Re[F(k, w)] relative to energy (uj) at k = Qxdw 
in Fig. 3(f), with the symmetric and anti-symmetric com- 
ponents of Re[-F] at k = Qxdw shown in the upper and 
lower panels, respectively. We find that F is predomi- 
nantly symmetric at Qp^w and primarily antisymmetric 
at Qsdw- the other hand, F appears to have compara- 
ble symmetric and anti-symmetric components at Qcdw- 
The antisymmetric F is consistent with SDW scenario for 
k II {tt,7t) [15,16], whereas for k || (tt, 0)/(0, tt), symmetric 
and anti-symmetric F may be attributed respectively to 
PDW [13,14] and CDW [17-21]. 

In Fig. 4(a) we illustrate the energy histograms of the 
SC and PC (or CO) features for H = 0, 2, 4.5 and 6 T. 
A strong spectral shift from SC a,t uj — Asc to PC at 
^ = Vco is seen with increasing iJ, together with the ap- 
pearance of a third subgap (SG) feature at a; = A' < Age. 
For comparison, schematic histograms for conventional 
type-II superconductors are shown in Fig. 4(b), which ex- 
hibit spectral shifts from an initial uj — Asc to a con- 
tinuous energy distribution lu < Asci with an additional 
peak appearing at w = in the H <C Hc2 and T Tc 
limit. The fraction of the spectral downshift is approx- 
imately given by {n(,'^Q/2)/{^/3a'^/4:), which is linear in 
H. The predictions in Fig. 4(b) apparently differ from 
our empirical findings in Fig. 4(a). In Fig. 4(c) we sum- 
marize the Gaussian fitting parameters to the histograms 
in Fig. 4(a). We find that the values of both Asc and Vco 
remain invariant with H, whereas the CO spectral weight 
and linewidth increase with increasing H. 

Discussion. — The occurrence of two energy scales 
Vco(> Asc) and A'(< Asc) inside vortices and the field- 
induced energy- independent wave- vectors Qpdw Qcdw 
and Qsdw strongly suggest that the vortex-state quasipar- 
ticle tunneling spectra in Y-123 cannot be explained by 
simple Bogoliubov quasiparticle scattering interferences 
alone [24-26]. On the other hand, these findings may 




0246 0246 0246 



H(T) H(T) H(T) 

Fig. 4: (color online) Field-dependent spectral evolution at T 
= 6 K in Y-123: (a) Energy histograms derived from STS 
data for H = 0, 2, 4.5, and 6 T, showing a spectral shift from 
Asc to Vco and A' with increasing H . (b) Schematic of the 
histograms for a conventional type-II superconductor in the 
limit of T <C Tc and H <C Hc2- (c) Gaussian fitting to the 
histograms in a reveals nearly field-independent Asc and Vco 
(top), decreasing SC and increasing CO spectral weight with 
increasing H (center), and increasing SC and CO linewidths 
with increasing H (bottom). 

be compared with the scenario of coexisting CO's and 
SC [13-21]. In particular, we find that the energy scale 
Vco and the wave- vector Qcdw rnanifested in the vortex- 
state spectra are consistent with the CO parameters de- 
rived from our Green's function analysis of the zero-field 
data. That is, we assume that the ground state of Y-123 
consists of coexisting SC and CO so that the corresponding 
mean-field Hamiltonian is given by the following expres- 
sion [19,20]: 

Hmf ~ Ti-sc + T^GO 

= $I^k4 „Ck,a - ^ Asc(k)(c[_^cLk,^ + c-k,iCk,T) 

k,a k 

+ Xl^C0(k)(c^Q^Ck,a +4,aCk+Q,a)- (2) 
k.Q 

Here Asc(k) — Asc (cos 20k) for d^2_j^2-wave pairing, k 
denotes the quasiparticle momentum, = ta.n^^ (ky/kx), 
^k is the normal-state eigenenergy relative to the Fermi 
level, and c are the particle creation and annihilation 
operators, and a —1,1 refers to the spin states. We in- 
corporate realistic bandstructures into ^k for direct com- 
parison with experiments [19,20]. By assuming a Fermi 
surface-nested CDW along (tt, 0)/(0, tt) as the relevant 
CO, we diagonalize TYmf and obtain the bare Green's func- 
tion Go (k, Lo) . The effect of quantum fluctuations may be 
further included by solving the Dyson's equation for the 
full Green's function G(k, a;) [19,20]. Thus, the quasi- 
particle DOS may be derived from G(k, ui) as detailed in 
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Refs. [19,20]. For finite temperatures, we employ the tem- 
perature Green's function [19,20]. 

Following the approach outlined above and detailed in 
Refs. [19,20], we can account for the zero-field tunneling 
spectra shown in Fig. 1(a) by the following fitting param- 
eters: Asc = (20 ± 1) meV, Vqo = (32 ± 1) meV, and 
Qco = (0.257r,0)/(0,0.257r) for the CDW. We note that 
the energy scale Vco = (32 ± 1) accounts for the shoulder- 
like features at w = ±Aeff ~ ±38 meV in the tunneling 
spectra, where Aes = [Al^ + V^^]'^/'^ [19,20]. Moreover, 
the magnitude of Vco = (32 ±1) agrees with that of the 
spin gap observed in neutron scattering data of optimally 
doped Y-123 [32] and the PG energy seen inside the vortex 
cores, as shown in Figs. 2(c)-(g). On the other hand, the 
corresponding wave-vector Qco is along the Cu-0 bond- 
ing direction, which remains in the vortex-state and is 
identified as the mode Qcdw the FT-LDOS. 

As an interesting comparison, we note that our recent 
vortex-state STS studies on an electron-type optimally 
doped cupratc Lao.iSro.9Cu02 (La-112) also revealed PG- 
like features inside vortices [34], except that the PG energy 
in La-112 is smaller than Asc- This finding may also be 
interpreted as a CO being revealed upon the suppression 
of SC. Furthermore, the smaller energy associated with 
the PG-like features inside vortices of La-112 is consistent 
with the absence of PG above Tc in electron-type cuprate 
superconductors [34] . These findings from the vortex-state 
quasiparticle spectra of La-112 are in contrast to those of 
Y-123 and the differences may be attributed to the differ- 
ent magnitude of Vco relative to Age [34]. 

Conclusion. — Our spatially resolved scanning tun- 
neling spectroscopic studies of Y-123 in the vortex state 

have revealed various novel spectral characteristics, in- 
cluding two energy scales (Vco and A') other than 
the SC gap Aco inside vortices and three accompany- 
ing energy-independent wave- vectors Qp^w? Qcdw and 
QsDW- These results cannot be reconciled with theories 
assuming a pure SC order in the ground state [35] . Rather, 
they are consistent with the CO scenario and suggest im- 
portant interplay between SC and various collective exci- 
tations in high-Tc superconductors. 
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